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I , Abstract 

Cn , We characterise the homogeneous and isotropic gauge invariant and 

quasifree states for free Dirac quantum fields on Robertson- Walker space- 
^~^ ' times in any even dimension. Using this characterisation, we construct 

,*li , adiabatic vacuum states of order n corresponding to some Cauchy sur- 

|i^ ' face. We then show that any two such states (of sufficiently high or- 

^— N I der) are locally quasi-equivalent. We propose a microlocal version of the 

VQ ' Hadamard condition for spinor fields on arbitrary spacetimes, which is 

("^ ' shown to entail the usual short distance behaviour of the twopoint func- 

^\ , tion. The polarisation set of these twopoint functions is determined from 

0^ ' the Dencker connection of the spinorial Klein-Gordon operator which we 

show to equal the (pull-back) of the spin connection. Finally it is demon- 
strated that adiabatic states of infinite order are Hadamard, and that 
those of order n correspond, in some sense, to a truncated Hadamard se- 
ries and will therefore allow for a point splitting renormalisation of the 
expected stress-energy tensor. 






X 



I Introduction 

In many cases of physical interest, for example the early stages of the universe or 
stellar collapse, one is naturally led to the problem of constructing quantum field 
theories on a non-static curved spacetime. Numerous papers have been devoted 
to the study of linear scalar fields on such backgrounds, but less has been done 
for fields with higher spin, mainly because the analysis of multicomponent fields 
is technically more involved. The aim of the present paper is to partly fill this 
gap for the case of a Dirac field on a curved spacetime. 

Quantum field theory in curved spacetime (in short, QFT in CST) is best 
described within the algebraic approach to quantum field theory, which started 
with the work of Haag and Kastler JJ^, for an overview see PHI- In this approach 
one deals with a net of C*-algebras {£^{ff)}ff<zM of observables localised in a 
spacetime region & C M . The algebra s^ = U^cAf ■S3?'(i^) is called the 'quasilocal 
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algebra'. Quantum states in the algebraic framework are positive normalised 
linear functionals on the quasilocal algebra. One of the major difficulties of 
QFT on CST is to pick out physically reasonable states. 

It has become widely accepted by now that, for a linear scalar field, the so- 
called "Hadamard states" are good candidates for physical states. These states 
are distinguished among other states by the particular form of the singular part 
of their twopoint function. A mathematically precise definition was given in |18j . 
The following facts about Hadamard states are known: 

They allow for a point-splitting renormalisation of the stress-energy tensor 
Tfj.1, |S3 ■ Verch [22] has shown that Hadamard states are locally quasi-equivalent 
and he has also shown local definiteness in the sense of Haag et. al. |12j . 
Radzikowski [27| discovered that (quasifree) Hadamard states, initially defined 
by the singular behaviour of the twopoint function in position space, can also 
be characterised by the so-called 'wave front set' of that twopoint function, a 
central concept in the mathematical subject called 'microlocal analysis', aimed 
at describing the singular behaviour of distributions. (Since these techniques 
do not belong to the daily used toolchest of the theoretical physicist so far, 
we give a brief introduction to this subject in the appendix.) The microlocal 
characterisation of Hadamard states is much easier to check in many cases where 
an explicit expression of the twopoint function cannot be obtained and has 
already led to important progress in the subject. It played an important role 
in the proof |2H| of Kay's conjecture in (axiomatic) QFT in CST, in the 
derivation of "quantum inequalities" ^ and for the perturbative construction 
of self-interacting quantum field theories in general globally hyperbolic curved 
spacetimes 0. 

There is a pre-existing notion of a Hadamard state for Dirac fields in a 
curved spacetime j^l OS] , analogous to the condition on the singular part of 
the twopoint function for a linear scalar field. It had been expected that there 
should also be a microlocal characterisation as for the scalar field, but the details 
had never been spelled out. One purpose of this paper is to close this gap. The 
microlocal condition that we propose is similar to that in the spin-0 case. But it 
differs in that, unless the state in question is assumed to be charge invariant, it 
needs to be imposed separately for the positive and negative frequency twopoint 
functions. We show ('Thm. llV.2|l that our microlocal notion of Hadamard states 
coincides with the concept based on the short-distance behaviour put forward 
in J19l \'^'^\ . Our result is the counterpart of similar theorem by Radzikowski |^7] 
obtained earlier for a scalar field. 

New questions also arise in the spin-1/2 case that have no counterpart in 
the spin-0 case. For example, it is natural to ask what the most singular com- 
ponents of the twopoint function are. It appears that the natural mathematical 
setting to analyse this question is provided by the concept of the 'polarisation 
set' 13 of a vector valued distribution (such as the twopoint function in the 
spin-1/2 case), a notion which refines that of the wave front set of a vector val- 
ued distribution. Making use of a theorem by Dencker ^ and the equations of 
motion, we determine the polarisation set of the twopoint function fThm. lTV.lll . 
corresponding to a Hadamard state of the Dirac field. Along the way, the prop- 
agation of singularities for the spinorial Klein-Gordon operator is obtained in 
Pron. lTvn 

On a Robertson- Walker spacetime and for free scalar fields, there exists the 



concept of 'adiabatic states', which was introduced a long time ago by Parker [^S] 
and put on a rigorous mathematical footing by Liiders and Roberts |55]. The 
main idea behind this concept is the following. If the scale factor R{t) in the 
Robertson- Walker metric is constant in time, then there is an unambiguous 
notion of positive frequency solutions to the Klein-Gordon (KG) equation, and 
one can use these to define a ground state. If R{t) is not constant, then no global 
ground state exists and the positive freqency solutions have to be determined 
dynamically off a given Cauchy surface (corresponding to the instant of time at 
which one wishes to define a vacuum- like state). They are usually found by a 
WKB-type ansatz and a subsequent iterative approximation process. Adiabatic 
states of order n (at the time in question) are then the states defined from the 
positive frequency solutions obtained after n approximation steps. Recently, 
Junker ^B] showed that the problem of finding the 'right' positive frequencies 
at an instant of time can be viewed as the problem to factorise the KG-operator 
into positive and negative frequency parts near the Cauchy surface in question. 

In this work we present a construction (cf. Def. IV.1|I of adiabatic states for 
Dirac fields on {N + l)-dimensional Robertson- Walker spacetimes, based on a 
factorisation of the spinorial KG-operator near a Cauchy surface. While the 
construction is similar to the scalar case, there are also important differences; 
for example, special care needs to be taken in order to obtain a manifestly 
positive state. In Prop. IV. II we explain how the factorisation of the spinorial 
Klein-Gordon operator (and hence the split into positive and negative frequency 
solutions) is related to the charge conjugation symmetry of the Dirac equation. 

In the scalar case, Liiders and Roberts have shown ^2 that the adiabatic 
states of sufficiently high order are all locally quasi-equivalent and Junker jl5| 
established that those of infinite order are of Hadamard type (which together 
implies that adiabatic states of sufficiently high order are locally quasiequivalent 
to a Hadamard state^). In Sec. IV of this work we show that these results also 
hold for our adiabatic states in the spin-1/2 case. Furthermore, we show that 
they correspond, in some sense, to a truncated Hadamard series and therefore 
allow for a point-splitting renormalisation of the stress tensor T^^ in the same 
way as scalar fields; in other words, our adiabatic states do not lead to infinite 
energy fluxes. 

Some of our results in the context of Robertson- Walker spacetimes can be 
generalised to arbitrary globally hyperbolic spacetimes, for example the con- 
struction of Hadamard states and a similar criterion for local quasiequi valence, 
based on the theory of pseudodifferential operators. For these and related issues 
we refer to a forthcoming paper. 

Concerning Thm. IIV.II A proof of this was given in an earlier version of the 
present paper which was, however, unfortunately incorrect. The first correct 
proof was given by K. Kratzert [2I|, see also |5J. These papers in turn biult on 
earlier unpublished work by Radzikowski "5^. We are grateful to K. Kratzert 
for communicating his results to us and for pointing out the error in an earlier 
version. In the light of his derivation of that result we were able to repair our 
earlier proof. It is included here since it is somewhat different from the proof 
in 1201 EH- 

^Originally, it had been claimed in ,jS that also adiabatic states of finite ordere were 
Hadamard. This has been corrected by the author of that paper in the meantime, cf. Iltil 



II The Dirac field on Robertson- Walker space- 
times 



The aim of this section is to recah the structure of the Dirac equation on 
Robertson- Walker spacetimes. We assume that the reader has some famiharity 
with the concept of spinors and the Dirac equation in curved space, as described 
e.g. in ref. |4j. 

The homogeneous and isotropic spacetimes in {N + 1) dimensions are of the 
form M'^ = E X E'^, the spatial section E"^ being the A^-dimensional sphere §^ 
for K = +1, the Euclidean space M^ for k = and the (real) hyperbolic space 
H^ for K = — 1. The line-element on these spacetimes is 



(1) 



where 



dsl = dt^ - R^{t)[dd^ + fl{e)dnlj_{\ 



{svaO for k = +1 
e for K = 

sinh 9 for k = — 1 

and dVLJ^_-y is the line-element on §^^^. The above spacetimes are models for a 
closed, flat or hyperbolic universe with positive, zero or negative curvature. We 
denote by n'^S^ = dt the future pointing unit vector field normal to the Cauchy 
surfaces E** and by /i^" — g'^" — n^^n" the induced (negative definite) metric on 

The spaces E" are homogeneous for the groups G+^ = Spin{N + 1), G" = 
Spin{N) X R^ and G"^ ^ Spin{N, 1) respectively, i.e. E" = C^/K, where 
K = Spin{N). We shall omit the superscript k when not necessary and assume 
that N is odd, iV > 3 in order to simplify the exposition. 

In order to bring out the 2 by 2 block matrix form of the Dirac equation in 
RW-spacetimes, it is useful to define the associated vector bundles (we view G 
as a K principal fibre bundle over E) 



E^ ^GXr< 



E^ = GXf^ 



,(JV-l)/2 



where t is the fundamental representation of K and f the conjugate of that 
representation. They are related to each other by ZT{k)Z~^ = f(fc) for all 
k € K^ where Z is some unitary matrix. The spinor and cospinor bundles 
(restricted to some Cauchy surface E(t) = E x {t}) are related to these bundles 
by DM \ E(i) = E^ ® E'' and D*M \ E(i) ^ E^ ® E^ . Decomposing 



C°°{M,DM) 3 iP 



X 



6,xeG°°(M,ii;^), 



(2) 



one can write the Dirac operator as the following 2 by 2 matrix operator: 



iN , 



{^ - m)^ = 7° ido + — 9o logi? 



where 



ogR- 


m 

iR-y 


iR-y 


/ 


X 


1 ■ 






-1 











The operator JJ is the Dirac operator on i?"^, defined without the scale factor R. 
The (generalised) eigenfunctions of this operator (normalised w.r.t. the natural 
inner product for sections in i?^), 'PXks ~ ^^^'^ks^ *"^^ ^"^ found in terms of 
special functions. The labels (fc, s) mean 



k = (fc, I, m) with 



J fc e N + N/2, 0<l <k- N/2, m = 0,--- ,di for k = +1, 



[ke 



1 = 1,2,..., m = 0. 



for K = 0,-1, 



and s = ±1. di is the degeneracy of the eigenvalue I + {N — l)/2 for the Dirac 
operator on S^~^, given e.g. in j31) . For an explicit representation of the 
functions Xks ^^ refer to 2, in the cases when k 
treated in the appendix. 



-1, +1. The case k = is 



In order to diagonalise the Dirac Hamiltonian, one defines the spinors 



ks 







uZ =R-^/^Uus 

ks 





^ks. 



where 



"'-71 



^/T+rn/uJk 



-S a/I - m/u) k 

/uk \/\ + m/wfe 



(3) 



is a unitary matrix and where Wfc = \Jvn?- + fc^/i?^ are the instantaneous fre- 
quencies of the mode. The spinor fields m- form a complete set of generalised 

ks 

eigenfunctions for the Hamiltonian 

and the helicity operator 
S = sign(is''V^) = sign 



m 



iR-^^^ 



iR~^j7 
—m 



if 



P — U|l/2gP''CT...p 



n^la 



■7p, 



(4) 



with energy iwfc and helicity s at the corresponding instant of time, 



ks 



ks ks 



ks 



(5) 



We next define the Dirac conjugate and the charge conjugate of a spinor or 
cospinor, which will be needed later on. Note that, as in any associated vector 
bundle, there is a one-to-one correspondence between sections x in E'^ ^nd 
smooth C^ -valued functions x^ on G such that x^igk) = T{k)^^x^{9) 

(a similar identfication can be made for sections in E"^). The charge conjugate 
(denoted by ip'^ or C'0) resp. Dirac conjugate ( denoted by ip or (3ip) of a spinor 
■0 with decomposition |5J are then defined by 



r = 



-(Z0^t)v 
(Zx^t)v 



^ = 



UAtW 



(3 is an antilinear map from DM to D*M and C is an antilinear map from 
DM to DM for which {C,H} = 0. The charge resp. Dirac conjugate of a 



cospinor is defined in a similar way. The action of the symmetry group G of the 
RW-spacetimes on spinors is defined as follows: Firstly, one has an action U on 
sections in E'^ , given by^ 

(f/(5)x)(^):=(x^(5-)r(^)- (6) 

That group action then extends, by means of the isomorphisms DM \ E(i) — 
E'^ S'^, to an action [/ = [/© t/ on spinors over M . An action of G on 
cospinors is defined by U{g) — (3U{g)(3~^. 

An operator B on C^ (M, DM) is called isotropic if it commutes with every 
U{g), g & G. By abuse of notation, we will also use this term for operators 
defined only on a single Cauchy surface. Such operators have a mode decompo- 
sition 

(A,s/.).=/dfc^5:6-:^4-^, (7) 

-' s pq 

where bks is some matrix valued function of the labels and /- = (/, ui )t, the 
scalar product on a Cauchy surface I](t) being defined by 

(/i,/2)t= / (/i7M/2)d5^ (8) 

The correspondence B ^^ b respects products and taking hermitian adjoints, if 
these are well-defined. For a proof of the above facts we refer the reader to ^] . 
For later use we mention that U can be viewed as a unitary representation of G 
on J^ = L^(E(t), Z5M), the space of square integrable spinor fields w.r.t. the 
above inner product. A bounded isotropic operator B on J^ corresponds to an 
essentially bounded function bks ■ 

The Dirac operator on a globally hyperbolic spacetime (such as RW-spacetimes) 
has unique retared and advanced fundamental solutions ^j^ and ^a, see 0], sat- 
isfying 

(ip - m)^A = (S'a(JF - m) = 1, {ip - m)^R = ^i?,(«F " "^) = 1. 

and (by J we mean the causal future resp. past of a region in spacetime) 

supp{^aI) C J+{supp{f)), suppi^jif) C J^{supp{f)) 

for any compactly supported testspinor /. p = ^a ~ ^r is called the causal 
propagator. 

Ill Local algebras for the Dirac field and invari- 
ant states 

III.l Local algebras of observables 

The Dirac field on globally hyperbolic manifolds can be quantised in a straight- 
forward manner. For convenience, we review the basic steps here, details can 

^ [/ is in fact the representation of G induced by the unitary representation r of the closed 
compact subgroup K. 



be found in 4 , which wc follow closely. As above, let Jift — i^(S(t), I?Af) and 
Jif-l its topological dual, identified with i^(E(t), D*M) (with the inner product 
defined in a similar way as in ^). The field algebra F is the uniquely defined 
unital C*-algebra CAR(J?^) generated by "time t"-field operators ^t{f) and 
^t{h), smeared with square integrable spinor fields / S Jit resp. cospinor fields 
h G J^l 1 which satisfy the "equal time i" anti-commutation relations (CAR's) 

All other anti-commutators are trivial. One also defines (iV -I- l)-smeared field 
operators by '!'(/) = '^t{^f \ S(i)), where / is now a compactly supported, 
smooth spinor field on M. In a similar way, one defines ^(ft), where h € 
C^{M, D*M). The field {N + l)-smeared field operators satisfy by definition 
the field equations 

^{{-iV - m)h) = ^{{ip - m)f) = 0, 

and the CAR's (see e.g. ref. |4)) 

{*(/), ^{h)} = i^{h, /)1 for aU / € C^{M, DM), h e C^{M, D*M). 

From this it follows at once that the definition of F is in fact independent of 
the choice of Cauchy surface made above. The algebras of fields localised in 
a spacetime region ^ are defined to be the C*-algebras F{0') generated by 
field operators smeared with test functions supported in ff. The algebras of 
observables localised in ff are given by £^{ff) = F{ff)™'^'^, where we mean 
the subalgebras generated by products of an even number of fields. From the 
support properties of the causal propagator, one can easily deduce that spacelike 
commutativity holds for the algebras of observables, 

K(^), i/(^')] = {0} if ^ and ff' spacelike. 

The group actions U and [/ of G on spinors resp. cospinors give rise, by standard 
results on the CAR, to an action by ^-automorphisms a^, g G G on the field 
algebra F. The action of these automorphisms on field operators is given by 

ag-^tU) = *t(f/(5)/), aa^t{h) = ^t{U{g)h). 

In the following, we will drop the subscript t at the "time i" -field operators. 
This should cause no confusion, as it will be clear what is meant from the 
context. 

III. 2 Invariant, quasifree states 

A state Ld on s^ — F'="«" is said to be isotropic \iLo{X) = Lo{agX) for all g G G 
and X e £^ . It is said to be gauge invariant and quasifree if there exists an 
operator < B < 1 on L'^ {T.{t) , D M) such that 

u:{^>{h) . . . ^{fn)^{h,) . . . *(/!„)) = 5„„det {{K ^/j)),,,=i,....„ : (9) 

h, e L^{i:{t),D*M), fj e i^(s(t), dm). 

The term "gauge invariant" refers to the fact that only monomials with the 
same number of ^ and ^ fields have a nonzero expectation value in the state 



uj. Clearly, a gauge invariant, quasifree state is isotropic if the corresponding 
operator B is. 

One can easily show that the GNS-construction {■Kuj,-^u,^u)) of a gauge 
invariant quasifree state gives the following: ^^^ is the antisymmetric Fock- 
space over .J^t ® ^', ri^ is the Fock- vacuum and the representation n^ is 



ir^{^{f)) = a+[B^'\^f)\nt)] + «-[(! - S)i/2(^/)r E(i)]*, / e C^{M,DM). 

Here, a± arc the destruction operators on .^^^ for particles and antiparticles 
corresponding to the respective copies of ^J^, which satisfy the usual anticom- 
mutation relations, 

{a+(/i)*,a+(/2)} - (/i,/2)t, {a-(/ii)*,a_(/i2)} = {hi,h2)t, 

(all other anti-commutators vanish) and a+{f)Qi^ = a^{h)fl^ = 0. 

We next want to ask when two given quasifree, gauge invariant isotropic 
states are locally quasiequivalent. Let ui and uj2 be two such states, corre- 
sponding to isotropic operators Bi and B2 (acting on some Cauchy surface 
I](t)) with decompositions 61 and 62 as in Eq. {Tj). 

Theorem III.l. The states lo\ and L02 are locally quasiequivalent provided 

ess.sup[(l -f |fc|)^+lfoifc, - &2fc.||] < 00 

(k,s) 

(we mean the matrix norm in C^J for some e > 0. 

Proof. Let us choose a region ^ of the form D{'Ta), where we mean the domain 
of dependence of some open subset 'rf with compact closure of the Cauchy 
surface S(i), i.e. the set of all x G J^C^) such that every past resp. future 
directed timelike or null curve starting at x hits '^. Let us first show that the 
restrictions of the states to a subalgebra £/{ff), G = D{'^) are quasiequivalent. 
Regions of this particular shape are convenient, because the algebras £^{ff) 
are then isomorphic to the algebras CKR{Jt^) constructed from the Hilbert 
space .J^ = L^{'!f,DM) which is a closed subspace of J^. The restricted 
states uJi^2 \ G then correspond to the operators E<gB\^2E<g on J^, where E'^ 
denotes the projection on this subspace. One then knows, by a well-known 
result of Powers and St0rmer [2S Thm. 5.1], that the states wi \ 6 and W2 \ & 
on the algebra CAR(J?^) are quasiequivalent if and only if 

\\(E^4'BxE^^f^'' - {E.^B2E^4'f'''\W^. < 00, 
\m^{I ~ Bi)E^)^/^ ~ {E^^il - B2)Ev)^/^n.s. < 00, (10) 



where we mean the Hilbert-Schmidt norm in J^. Using the Powers-St0rmer 
inequality 



\A,^A2\\l,s.<\\Al-Af 



valid for any two bounded operators (we mean the trace norm), one concludes 
that Eqs. ^^ hold provided 

WEvBiE^ - E<^B2E^\\tr < 00. (11) 



As above, let H be the Hamilton operator on J?^. Trivially, we can write {p E R) 

E^B.E^^ - E^B^E^ = {E^\H\-P'^)[\HY>'\B^ - B2)\HY>'WH\-^'^E^). 

By assumption, the operator \HY'/'^{Bi - i?2)|-ff|^/^ is bounded ior p < N + e, 
therefore Eq. l(TT)) will hold if E<g\H\~'''/'^ can be shown to be Hilbert-Schmidt 
for such a p, since the product of two Hilbert-Schmidt operators is in the trace 
class. To see this, let us pick an orthonormal basis {/n}neN of spinors in J^. 
Then 

n 

= Y.j'^^ll(^"iR"+^"r""\fl-J- (12) 

n qs 

In order to estimate the r.h.s. of Eq. (|12|l . we exchange the dfc-integration and 
the summation over n (this is justified, because the resulting expression turns 
out to be absolutely convergent). We have. 






The sum over g, s, fc at fixed fc of the integrand is independent of x and equal 
to twice the spectral function Pj^ik), defined by 

■FW(fc) = y^Xfcim^(0)^Xfcims(0). 
Ira 

therefore we have found 

||ii;^|ifrP/2|l?i.g. = 2volC^)y(fcVi?' + m2)-P/2F^(fc)dfc. 
The spectral function Pjv is given by the following expressions 

'^{N-\Ml {N+k-\)\ for K - +1 

^ fe!(7V-l)! torK--|-i, 



I fc™~^ 

Pn{^) = < 2"/2vol(S"-i)r(Ar/2)2 

r(Af/2+ifc) 
r(Af/2)r(l/2+iA:) 



-^w- 



2 



for K = 0, 
for K = — 1. 



For K = — 1,+1 a derivation of these may be found in [2], the expression for 
K = is derived in the Appendix. P;v(fc) grows as k^~^ for large k for for all 
the homogeneous spaces E = R^,S^,]HI^, ensuring that £'<g'|-ff|~^/^ is Hilbert- 
Schmidt for p = iV -I- e. We have therefore shown that uji [" i^ is quasiequivalent 
to 0^2 t ^ for a-ny set ^ of the form D{^). Since it is enough to verify local 
quasiequivalence on a cofinal set of open subsets (such as the set of regions of 
the type D{'^)), this then proves the theorem. D 

IV General properties of Hadamard states for 
the Dirac field 

The plan of this section is as follows: We first give a microlocal definition 
of Hadamard states for Dirac quantum fields. After that we determine the 



polarisation set of the twopoint function of such states. FinaUy, we explain how 
our microlocal notion of Hadamard states is related to pre-existing notions of 
Hadamard states for the Dirac fields based on the short-distance behaviour of 
the twopoint functions. The definitions and results in this section apply to the 
case of a general, globally hyperbolic spacetime. The reader not familiar with 
the technical ingredients of the definition will find some notation and results 
from microlocal analysis in the appendix. 

The spatio-temporal twopoint functions of a state w are denoted by 

(p^+\K I) ■■= uji^imm, (p^-\K I) ■■= ^(^(/i)*(/)), 

where / € Cg°{M,DM) and h e C^{M,D*M). They are assumed to be 
distributions. We introduce the following (standard) notation: elements in T*M 
are denoted by (x, ^). We write (xi, ^i) ~ (a;2, ^2), if xi and X2 can be joined by 
a null-geodesic c such that ^1 — c(0) and ^2 = c(l), or if (a;i,^i) = (a;2,^2) and 
^1 null. We write xi >- X2 resp. xi -< X2 if the point xi comes after or before 
X2 according to the parameter on this curve. Moreover, we shall write ^ i> if ^ 
is future-directed and ^ < if it is past-directed. 

Definition IV. 1. A quasifree state uj for the Dirac field is said to be 'Hadamard' 
if 

WF'(^(±)) = {(xi,ei,X2,6) e r*M\{0} X T*Af\{0} | 

(xi,Ci)~(x2,6), Ci>WO}-:C(±) (13) 

Remark IV. 1. It follows from what we say in the proof of Thm. IVI.ll that 
the apparently weaker condition WF'(G^^^) C C*^^-* already implies equality of 
these sets. 

The conditions on the wave front set of ^^"''^ and ^^~' are in general inde- 
pendent, i.e., there exist states for which ^(+^ has the desired wave front set 
but for which fi^ ' has not. However, it is easy to see that this cannot happen 
for states u which are invariant under charge conjugation. By this we mean 
that Lo o ac = u}, where Uc is the *-automorphism defined by 

Below (Thm. IIV.2(I we prove that the above condition on the wave front sets of 
both twopoint functions implies that their singular behaviour in position space 
is that of a Hadamard fundamental solution. This result will in general not hold 
if only one twopoint function has the required wave front set. 

We now analyse in more detail the microlocal singularity structure of the 
twopoint functions by applying to them the propagation of singularities theorem 
Thm. IVIII.ll taking P in that theorem to be the spinor Klein-Gordon (KG) 
operator, given by 



1 

4^ 



P = n + -i?-t-m^ 



where R is the curvature scalar. This theorem is applicable because 

(P®1)^(=^) = (l(g)P*)^'=^) =0, (14) 

10 



which is in turn a direct consequence of the Lichncrowicz identity 

P ^ [ip - m){-i'p - m) 

and the fact that ^(^' satisfy the Dirac equation. By * we mean the transpose 
of an operator acting in the dual bundle (the bundle D*M in the case at hand), 
obtained from the natural pairing between spinor and cospinor fields on M . 
(The assumption in Thm. I VIII. II that P be of principal type (cf. Def. IVIII.3P 
is fulfilled, because P has a metric principal symbol po(a;,0 = —g^'^{x)£,f_t£,iy) 
In order to apply Thm. IVIII.ll we must calculate the Denckcr connection Dp 
associated with P, defined in Eq. (I39|l . In the case at hand, Dp naturally acts 
in the vector bundle n*{DM (g) Q^/^) \ Qp, where fl^/^ is the line bundle over 
M of half-densities, n : T*M -^ M is the projection map and 

Qp - {{x,0 I g'"'{x)U- = 0} C T*M. 

Proposition IV. 1. The Dencker connection Dp for the operator P is the par- 
tial connection in the pull-back of the vector bundle DM x J7^' ^ to Qp C T*M 
given by 

Dp^ix^^oTr*{W + -d\og\g\). 

Here, Xp^ is the Hamiltonian vector field over T* M corresponding to pq(x,^), 
jdlog\g\ is the natural connection in H,^'^ and ix^ is the insertion operator. 

Proof. The Dencker connection can be calculated from Eq. 1)39(1 , taking po — I 
(this means that q = po in that formula). Introducing an orthonormal frame 
e^da;^ for the metric tensor, rj""^ = Q'^'^^u^t^ ^^'^ ^^^ going to local coordinates, 
we find 



{p„,Po} = 0, p%x,0 = 2^S^{x)^^ - zel{x)d,e^,{x)C 



and 



d_ 



Xp„(x,C) = -2e^^+2eUx)d,e^.{x)Ci 



where S^ = |[7",7'']eQV^ey6 is the spinor connection. According to Eq. l|S^. 
this means that 

Dp - -2^^^^ + 2e^(a;)a^e,.(a;) - 2S,{x)e + e;ix)d,e^ix)C ■ 

That this expression agrees with the geometric expression for Dp in the state- 
ment of this proposition now immediately follows, because V^ — 9^^ -I- S"^, 
|;9^1og|5| = — ie°i9^ej^, and because the differential of the projection map 
is given by dT:{-^) = ■^, d-K{-^) = 0. D 

Having derived this tool, we can establish the following result on the polar- 
isation set of the twopoint functions of a Hadamard state. 

Theorem IV. 1. Any Hadamard state has the following polarisation set: 

WF;„;(^(±)) = {(xi,a,a:2,6,w^) I (2^1,^1,2:2,6) e C^^\w^ D,,M(g>D^,M, 
and w c'^{xi,X2)b = ^ii g for some X £ C}. 
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Here, unprimed spinor indices refer to the point xi, whereas primed ones cor- 
respond to X2 and J^(a;i,a;2) is the bispinor of parallel transport in the bundle 
DM along a null geodesic joining xi and X2. The sets C^ ' had been defined in 

Eq. (frvlTi . 

Proof. We aim at using the propagation of singularities theorem, Thru. IVIII.ll 
combined with a deformation argument due to FuUing, Narcowich and Wald [Jj , 
first apphed in a similar context in |10) . 

By Thm. IVilL II and Eq. i|14|l . the polarisation set of ^'■^^ must be a union 
of Hamilton orbits corresponding to the operators P ® 1 and 1 (g) P* . By the 
Prop. Hvm sections over Qp, annihilated by Dp are pull-backs to T*M of sec- 
tions in DM over nuU-geodesics which are parallel with respect to V. Therefore, 
two elements {xi,^i,X2,£,2,w) and (a;i,i^i,a;2,^2j''"') of 7r*(Z?Af Kl I?*M) are in 
the same Hamiltonian orbit if 

for some A G C. Now let x € M and S be a Cauchy surface of M through 
that point. Then there is there is a convex normal U oi x and a convex normal 
neighbourhood V of E, containing C/, such that there is another spacetime 
(M, g) with Cauchy surface E and a corresponding causal normal neighbourhood 
V with the properties that: (a) {V,g) is isometric to {V,g) and (b) M contains 
a Cauchy surface Si and a a convex flat neighbourhood C/i contained in a 
convex normal neighbourhood Vi of Si such that D{Ui) D U (we mean the 
domain of dependence), where U corresponds to U under the isometry By the 
propagation of singularities theorem, it will be enough to show that jf?^^^ \ UxU 
has the desired polarisation set, because any pair of null related points can be 

transported along a null geodesic into a region of that kind. Let p \ V xV 
be the pull-back of the twopoint functions to the deformed spacetime {V,g). 
By the propagation of singularities theorem and the equations of motion on 
the deformed spacetime it will induce a Hadamard distribution on all of M . 

Furthermore (p \ {j x U will have the required polarisation set if fi \ 
Ui X Ui has, again by the propagation of singularities theorem. But Ui C M 
is contained in a flat portion of spacetime, so effectively our theorem has to be 
shown for Minkowski space only. 

So let ^J„j„]j be twopoint functions of a Hadamard state in Minkowski space. 
By our Thm. IIV.2I all Hadamard states differ by a smooth piece only, so we 
might restrict attention to the vacuum in Minkowski space, 

where A^j^jj^j, are the ordinary positive rcsp. negative frequency twopoint func- 
tions for the KG-operator r/'^'^d^di, + m? in flat space. It is not difficult to see 
from the definition of the polarisation set and the definition of ^„jj„i^ that one 
must have 

WFpoK^iuL) C {(xi,C,a;2,-e,w) I (a:i,^,X2,-0 e WF(A^^^,), 

w^cl+piov some c £ C and /3 € C^+^j. (15) 
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Now 



M -^)® i]^SL = ^LILii ® {-^P - '^)\ = 0- 



therefore, using that ^ is a principal symbol of i^ + 771, one can conclude (by the 
definition of the polarisation set) that 

lw = wi = 0, 

where (xi, ^, xi, — ^) e "^'^ {p'ai\a\i) ■ Since the form of w is already restricted by 
Eq. H15|l . it is easy to see that these equations imply rf-'^i,p.iiv — and c = 0. 
From the equation ^j^ = it follows 

tr(7"^ . . 7"«-7^+2e„,...„„_,.,/?^) = 

where £^...2^ is the totally antisymmetric tensor in A^ + 1 dimensions and 

Using standard identities for traces of gamma matrices and 

we find that ^p/3o- — (,aPp = 0. Since ^ y^ 0, this implies p — w = \- ^, which in 
Minkowski space is just the condition on the polarisation that was claimed. D 

In PSlEj, the authors give a definition of the Hadamard condition for Dirac 
fields in terms of the singular behaviour of the associated twopoint functions in 
position space. We now investigate the relation between the two definitions. A 
local version of the definition of |23| EI may be stated as follows. Let ^ be a 
convex normal neighbourhood in AI . On ff, one defines the bidistributions 

(N-l)/2 „+i 

where a is the signed squared geodesic distance between the points xi, X2 G ^, 

ae{xi,X2) = ct{xi,X2) + 2ie{t{x2) — t{xi)) + e^, e > 

and t is some global time function. As usual, we mean the bidistribution ob- 
tained by smearing with smooth spinor fields first and then taking e to zero. The 
bispinors Uj , Vj are determined recursively by the N + 1 dimensional analogue 
of the Hadamard transport equations [23 for the spinorial KG-operator P and 
depend only the geometry of the spacetime in ff. By construction, 

(P (g) l)i?i=^) = {1(E) P*)i?^=^' = mod C"* 

on ^. The local version of the global Hadamard definition'^ given in [SSI E| 
is as follows: A state is said to be locally Hadamard if its associated twopoint 
function satisfies 

^(±) = ±{-iy7 - m)H^^^ modulo C" on ff for all 77. (16) 



^Not all points in M may be connected by a unique geodesic line, therefore the definition 
needs to be refined if one also wants to exclude spacelike singularities. We refer to |18| for 
details. 
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In fact, it is only necessary to require this equation for '+' (as the authors do 
in ref. ISSLIlOp . since it automatically follows from this that 

(£•(-) = i{-if - m)E - ^'^+^ 

= {-ip - m) {iE - iJi+) ) mod C" 
= -{-if-m)Hl-^ modC", 

where we have used that 

i7,(+) - H^-^ = iE mod C"+\ 

and where E is the causal propagator for P. This follows from the fact [HI that 
E contains the same coefficients Uj , Vj as Hn , multiplied only by different 
singular pieces, which, as it is well-known in QFT, may be combined to give the 
above identity between the different types of propagators. 

We now show that Eq. (|16() follows from our microlocal definition. This 
observation has been made first by Radzikowski "57" for the scalar field and his 
proof can be adjusted to the spinor case as well, although the situation is more 
complicated. 

Theorem IV. 2. Let ui be a Hadamard state in the microlocal sense of Def \iV.l\ 
Then its twopoint functions also satisfy Eq. H16|l. 

Remark IV. 2. It can also be shown that the global version of the Hadamard 
condition IjlGfl (as for example spelled out in |33[ I19j) implies the microlocal 
Hadamard condition. A proof of this is almost identical to the proof of the 
corresponding statement for the scalar case. We refer to refs. |19II2()| for details. 

Proof. Let us denote by Ea, Er, Ep, Ep the advanced, retarded, Feynman, anti- 
Feynman parametrices of the spinorial Klein-Gordon operator P. They are 
known to be determined (modulo C°°) by the equations 

PEj = 1 mod C^, j ^A,R,F,F 

and their wave front sets |13l Thm. 6.5.3]. For the Feynman and anti- Feynman 
parametrices, these wave front sets read: 

WF'(i?p) = {(xi, a, 2^2, 6) I (2^1,^1)^(2^2, 6), 

a i> if xi >- X2, a < if xi ^ X2 } (17) 
and 

WF'(i?^)-{(a;i,ei,X2,6)|(xi,Ci)~(x2,6), 

a > if xi ^ X2, a <i if xi y X2 }■ (18) 

Actually, in ^J Thm. 6.5.3] only the case of a scalar operator with metric 
principal part is treated. Inspection of the proof however shows that it may 
be extended to operators with metric principal part acting in vector bundles 
such as P. We also need the advanced, retarded, Feynman and anti-Feynman 
parametrices for the Dirac operator, given by 

^, = i-if-m)E„ j^A,R,F,F. 
14 
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By the anticommutation relations, one infers that 

^(+) + ^(-) = ^^ = ^{^^ - ^j,). 

Let us define 

^F := 10+^ + ^A = -i^^-^ - 

fp := -*^(+) + ^K - *^(-^ + ,^A. (19) 

Our aim is now to prove that 

modulo a smooth kernel. To this end, we first show that 

WF'(^f) C WF'{Ef), WF'(^^) C WF'{Ep). (20) 

In order to see why this must be true, consider a point (xi, ^1,2:2,^2) in wave 
front set of fjp such that xi ^ J~{x2)- Then, because ^a must be zero for 
such points by the support properties of Ea, it must hold that (xi, ^1, X2, C2) G 
WF'(^(+)). Since (by the microlocal Hadamard condition, Eq. ITSJ) WF'(^(+)) C 
C'+' we find that (xi, ^1, X2, C2) can be in the wave front set of f!p if and only 
if ^1 1>0. A similar reasoning can be applied for xi e J~ {X2), this time using the 
representation (pp = —ip^' -\- ^r and exploiting the microlocal Hadamard con- 
dition, Eq. (|1H|I , for p^^\ Altogether one concludes from this that (xi , ^i , a;2 , ^^2) 
is in the wave front set of ^f if and only if (a;i,^i) ~ (^2,^2) and ^1 > for 
xi >~ X2 resp. ^1 oO if xi -< X2, which is just the set Eq. l(T7|) . We have therefore 
shown the first inclusion in Eq. I|20() . The second inclusion is treated in just the 
same way. 

Now, by definition, we have flp + pp = ^a + ^r- Applying the operator 
—ip — m to the relation 13, II, Eq. 6.6.1] 

Ef + Ep=Ea + Eii modC"", 

we find from this that fJp — ^f = Pp — ^p modulo smooth and hence that 

WF'(^f-^f)-WF'(^^-^^). (21) 

We had aheady shown that WF'(jS'f) C WF'(^f) and we also have WF'(^f) = 
WF'((— i^ — m)Ep) C WF'(i?F) (since the wave front set of a distribution 
cannot become larger when acting upon it with a differential operator), so the 
set on the left had side of this equation is contained in the set WF' (i^F)- By 
the same arguments the set on the right hand side Eq. (|21|l must be contained 
in W¥'{Ep). It therefore trivially follows that 

WF'(^F - Pf) c As, , ^Y'ipp - ^p) c Aqp , 

where 

Aqp = WF'(ii;F) n WF'(ii;^) = {{x, ^ x, o I ^^''(a^)^^^. = o}. 

Now by definition 

iP®l){pp-^F) = (1®P*)(^f-^f) = 
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modulo smooth. We are thus in a position to use the propagation of singularities 
theorem and we conclude that polarisation sets of the distribution (pp — ^p (and 
hence the wave front set) must be a union of Hamiltonian orbits for the operator 
P. Using the same arguments as in the proof of the preceding theorem, it is 
then easy to show that if (a;,^,a;,^) is in WF'(^i? — ^p), then this set must 
also contain nonzero vectors away from the diagonal, a contradiction. Hence 
^MY'{pp-pp) =0, i.e., 

pp-ppe C°°, 

as we wanted to show. Inserting this into Eq. (|19|l one gets 

i^(+) = pp-pA = {-if-m){Ep~EA) mod C°°. (22) 

It can be extracted from the analysis of the propagators in [HI El that 

Ep-Ea = iH^^^ modC"+^ 

This holds because all propagators have the same structure as Hn , i.e. the 
same functions Uj, Vj multiply different singular parts. These can be combined 
to give the above equation. Combining this equation with Eq. H22|l then proves 
the theorem. D 



V Adiabatic states 

V.l Definition of adiabatic states 

For the remainder of this work, we restrict attention to Dirac fields over RW- 
spacetimes. It will be clearly indicated if a result has a wider range of validity. 

We begin by constructing adiabatic states for the Dirac fields on RW-spacetimes. 
By term "adiabatic" we mean that any of these states should give a reasonable 
mathematical description of the concept of 'empty space' in the very small, 
i.e. in spacetime regions which are very small compared to the curvature ra- 
dius. The main ingredient in the construction is a factorisation of the spinorial 
Klein-Gordon operator into positive and negative frequency parts, which has 
been considered before in ^^, for the spin-0 case. (Such a factorisation is pos- 
sible on every globally hyperbolic spacetime with a compact Cauchy surface.) 
Let us outline our construction. Suppose one has constructed a 'Gaufiian fo- 
liation' / 9 t I— > S(i) of a neighbourhood of a globally hyperbolic spacetime 
M by smooth Cauchy surfaces S(t), / = [ioj^i] denoting some time inter- 
val. By 'Gaufiian', we mean that the vector field dt = n^dfj. associated to this 
foliation is geodesic and orthonormal to the Cauchy surfaces, in other words 
Vj^n^ — V^rij/ = and n'^'Vfj,n'^ = 0. It is then possible to make the decompo- 
sition (P is the spinorial KG-operator) 

P = - (m^ V^ + iK + T) {mf^Vf, - T) mod 0P"°° , (23) 

where T is a pseudodifferential operator (PDO) in OP^(S x /, DM) with prin- 
cipal symbol 



ai{T){x,0 = ^-h^''{x)U^^ (24) 
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acting 'surface- wise', i.e. arises as a smooth family {T{t)}tei of operators 
acting on the surfaces S(i). Here K = '^ f^n^ is the extrinsic curvature and 
H{t) — —ifih^'^jfj^Vi, + fim is the Dirac Hamiltonian. (For the various classes 
of operators and symbols, we refer the reader to the appendix.) Solutions T 
to Eq. (|23(l can be found by inserting the asymptotic expansion for the symbol 
of the sought-for operator in that equation and then determining the terms in 
this expansion iteratively. The iterative procedure may be stopped after a finite 
number n of steps, yielding operators T„ which differ from T by an operator of 
class OP"". In the case of a RW-spacetime, these will be isotropic. A more 
detailed discussion of how this works on such a spacetime and how the operators 
T and T„ are defined in that special case is given below, since we do not want 
to interrupt the present line of argument. 

Now fix a t G [toi ^i] and define 

L±(t) = r(t)±iJ(i), 

and similarly operators Ln,±{t) by taking T„(i) in this formula. Our definition 
of adiabatic states is based on the following lemma, valid on RW-spacetimes. 

Lemma V.l. Let n e N U {oo}. Provided R{t) ^ 0, then the operators L±{t) 
can be modified, smoothly in t, by a PDO with smooth kernel such that there 
exist hermitian, isotropic, positive operators Q{t) G OP^" on L^{'S{t), DM) 
satisfying (we omit the reference to the time t) 

L+QL*^+L*_QL^ = 1, (25) 

In the same way, if n is a natural number there exists Q„, positive and isotropic, 
satisfying this equation for Ln_±. 

Note: If M is a spacetime foliated by Cauchy surfaces, then the hermitian 
adjoint of an operator acting 'surf ace- wise' (such as in the lemma) is defined 
w.r.t. to the inner product Eq. ^ on each surface. 

Before we prove the lemma, let us give the definition of adiabatic states for 
the free Dirac field on a RW-spacetime. Let us set 

Bn — Ln^ + QnLn^+- (26) 

Then i?„ e OP°(S(i), Z5M) is hermitian, isotropic and by the lemma fulfills 
< S„ < 1. 

Definition V.l. The gauge invariant, quasifree states cj„ defined by the oper- 
ators Bn on L^(I](i), DM), n e NU {00} are called adiabatic states of order 
n at time t. 

It is straightforward from the definitions that the twopoint functions of such 
states are given by (suppressing the subscript n for the moment) 



P^+\h,f) = {{tn'^V,+T{t))fh,Q{t){in^^V^ + T{t))pf)^, 
P^-HhJ) = ((m^V^-r(i)*)^/j,Q(i)(zn^V^-r(t)*)^/)^. (27) 
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Construction of the operators T and T„ in the spin-1/2 case in RW- 
spacetimes: We have, 

P = -{if + m){iJ7 -m) 

= -{ini'V^ + iK + H){in>'V^-H), 

where we have used that in^^V ^ — H = fi{ip — m) and the identities 



w 



W^fi = 0, /i'"'7pV,^ = K, 



which foUow easily from n^Vpn"^ = 0, Wf^n^ — Vun^j. = and Vp7i/ = 0. 
Therefore Eq. H23|l is equivalent to the operator equation 

H^ + iKH+[in^'Vi„H]=T^ +iKT+[in^'Vi„T] mod 0P~°° (28) 



for all t G I, where T is required to have a principal symbol y—h^^^^^- While 
the above said holds for any globally hyperbolic spacetime with a Gaufiian foli- 
ation, we now specialise the discussion to RW-spacetimes, where one of course 
takes the obvious foliation by the homogeneous surfaces S'*. We consider sym- 
bols J X M 9 (i, fc) H^ bksit), taking values in the complex 2 by 2 matrices and 
carrying an additional helicity index s = ±. We will often omit reference to the 
matrix resp. helicity indices and simply write & € S" when we mean a matrix 
valued symbol of class S"(/,1R) (S) M2(C). Any such symbol b defines and oper- 
ator B via Eq. (0). Furthermore, if fo G S" then B G OP"(S x /, DM). To see 
this, one can argue in just the same way as in '15' , where a similar statement is 
proven. 

We wish to write T in the form Eq. Q with some matrix valued symbol r^s 
with principal part equal to uJk- Inserting the ansatz Eq. Q into Eq. (|28|l . one 
obtains the following 2 by 2 matrix system of ordinary differential equations 
with parameter k, 

^r + i^^T -f [r, d\ + t" '"° ="°° ih + q^h + [h, d] + e -: r, (29) 

where hks = diag[a;fc, ~ijJk\- To arrive at this equation, we used Eq. jSJ and 



where 



dZ = 



r\ JskmB, 

JskmR r\ 

2{ni^K^ + k-^) ^ 



(30) 



The drop-off properties in k required from the the matrix t^s uniformly in the 
time interval /, imply that the initial data must be carefully adjusted, and 
cannot be freely chosen. 

We try to find r as an asymptotic expansion of its symbol (e.g. in a certain 
sense 'in powers of A:'). First note, that since Eq. I|23|l should hold only up to 
the addition of a arbitrary operator of class OP^°°(S x /, DM), t need only be 
defined up to S^°°. It is therefore enough to find an asymptotic expansion for 

j 
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where "dj £ S ^•' and matrix valued. In order to get an operator with the 
right principal symbol we set ^o,s/c(0 = '^fc(i) for all t G I. We then define dj 
successively in such a way that the n'th partial sum r„ = ^" "dj in Eq. 131(1 
solves Eq. l(^ modulo S^~". This can be achieved by putting 

^„+i = ^ (*r„ + [r„, d] + r,^ + ^t„ - r) E S-"(/, R). (31) 

Following this procedure we can calculate terms of arbitrary high order in the 
asymptotic expansion for r. Any symbol with this expansion will give rise to an 
operator T factorising the spinorial KG-operator modulo 0P~°°. Moreover, the 
partial sum r„ obtained after n iterations will give rise to an isotropic operator 
Tn, which will solve Eq. H23|) up to OP^". T as well as r„ have the same 
principal symbol as \H\, i.e. Eq. H24I) holds true. 

Clearly, if the existence of operators Qn as in the lemma is known for a 
general globally hyperbolic spacetime, then one can still define operators -B„ by 
Eq. 1)26(1 . and these in turn give quasifree, gauge invariant states aj„. 

Proposition V.l. Suppose {M,g) is a spacetime with a compact Cauchy sur- 
face S, and suppose there exist operators Q„ fri g N U {oo}^ as in the above 
lemma. Then there are operators < -B„ < 1, differing from Bn only by an 
operator of class OP^"^ , such that CBnC~^ = 1 — i?„ for all n. 

Remark V.l. It follows from the definition of the charge conjugation auto- 
morphism that the states a)„ corresponding to i3„ are charge invariant. 

Proof. We treat the case of infinite order first and suppress the subscript n. We 
also introduce the notation A'^ = CAC~^ for operators A acting on spinors. We 
need to show that B can be modified modulo C°° to an operator B such that 
B'' = 1 - B and < S < 1 . Taking the transpose of Eq. ^^ we obtain 

P* = -(-m^V^ - iK - T^)i-in''V^, + T*) mod 0P"°°. 



It is not difficult to see that Pf = P*/ and T*/ = T*f. From this it follows 
that 

P= -{in''V^ + iK-T*){inf'V^ + T*) mod 0P"°° (32) 

But multiplying Eq. 1(23(1 with C from both sides and using that (m'^V^)'^ = 
—in^Wfj., we also have that 

P = -(in^V^ + iK - T''){in''Vf, + T") mod 0P"°°. 

Therefore, since the principal symbols of T* and T'^ are equal and since the 
factorisation is unique modulo 0P^°° once this information is known, we have 
shown that 

T*^T' modOP"°°. 

We redefine T by i(T + r*'=) and Q by liQ + Q") and L± by inserting the mod- 
ified definition of T. These redefined operators will then satisfy (remembering 
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that that H'^ = —H and using that A" = A*'^ for any operator A actmg on 
spinor fields over S(i)'*.) 



Q' = Q, Q > 0, L^ = LX 



and 



where A G OP °°. Since the left hand side of this equation is invariant under 
charge and herniitian conjugation and positive, we find A'' — A, A* — A and 
1 + j4 > 0. Since I](i) is compact, A is a compact operator and the projectors 
on all nonzero eigenspaces have a smooth kernel. Let F be the projector on the 
(finite dimensional) kernel of 1 + A. Then R = l + A + Fis strictly positive, 
[F, R] ^0, FR = F, F" = F,R''^ R. Let us write 

L± = R-^I^L^R-^'"^, Q = R^/^QR^/^. 

Then L^. ~ L*, Q >0, and the operator B defined by 



B = L+QL* 



^F 



satisfies B ^ B modulo smooth, B'^ ~ 1 — B and < B < 1, providing us thus 
with a modified operator with the desired properties. For arbitrary n > I, one 



proceeds in a similar way, this time using that T^ = T^ modulo OP 



D 



Proof of lemma IVTI 

Proof. The argument establishing the existence of Qn as in the lemma is the 
same for all n > 1, therefore, to lighten the notation we will only treat the case 
n = oo and drop the reference to n. We set £± = TzLh. The £± are then related 
to the PDO's L± in the statement of the lemma by Eq. ^. One finds after the 
first iteration, 



C+.fcs 



- 







Wfc 







modulo S^ . Further iterations change r only by symbols of order less or equal 
— 1 and will therefore not affect the above form of £± . Only the above form of 
the £± is used to argue the existence of Q, therefore the adiabatic order n is 
not important for our argument, as long as n > 1 . The proof of the lemma 
then amounts to show that t can be modified by a symbol of class S~°° such 
that one can find a, q G S~^(or rather one for each helicity), taking values in the 
complex 2 by 2 matrices, such that (* denotes the hermitian adjoint of a matrix 
and 1 the identity matrix) 



^+,fcs9fes^+,fcs + ^'^,ksQks£-,ks — 1, qks 



Iks 



Qks > Vfc > 0, S = ±. (33) 



The PDO Q corresponding to q by Eq. then obviously fulfills the claim of the 
lemma (note that the integral/sum in Eq. (0 is over positive k only). Taking 



^Here one must use that the foliation is Gaufiian. 
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the matrix adjoint of Eq. (|33|) . one observes that q can be taken to be hermitian. 
One may regard Eq. 1)331) as a hnear equation for q at each value of k and the 
hehcity index s, and thus write it as a 4 by 4 matrix system for the matrix 
entries of q 



(Iks 




'1 


22 
(Iks 




1 


ill 







fkl. 








Mks :?! - n ' (34) 



where Mks is a 4 by 4 matrix determined from the entries of i±^ks- Only using 
the above from of i± and the fact that S^S" C S"+™ one finds from Eq. (^ 



M = D{1 + A) where 
Dks = 'idiag[u;lu;l^^R-^/^^tiR^/'cJk),^R^''/^^tiR^/^cJk)], AeS 

Hence, if R{t) ^ 0, then Mks has a matrix inverse in S^^ for large k, 



-1 



r-l _ n-1 '^^ 

ks /_ J 

m=0 



^iks^^ksl^ ^ks 



Eq. 1)341) may therefore be inverted for large k and gives us a solution q to 

Eq. ()33)l . It follows directly from the above form of M~^ that qks has diag[(2a;fc)~^, i^LiJk)' 

as a principal symbol, therefore qks is positive definite for large k. We have thus 

constructed a solution to Eq. ()33)l if k is greater than some /co > 0. To find such 

a q also for < fc < fco: we may redefine 

^.^ ^, fT(i,/co) forO<fc<fco, 
1 T(i, /c) for k > ko, 

and arbitrarily for fc < 0, since T by definition does not depend on r(i, fc) for 
fc < 0. By what we have already shown, such a r trivially allows for a hermitian, 
positive solution of Eq. (|^ . but it is not yet a symbol (because its dependence 
on k is not smooth). We might however change the above definition of r in an 
arbitrary small neighbourhood of fcg to make it smooth (and hence a symbol), 
without making the corresponding matrix M singular. As we mentioned earlier, 
the resulting matrix q is automatically hermitian. By easy arguments based on 
the continuity of the construction, it will also remain positive, if that change is 
made arbitrarily small. D 

VI Properties of adiabatic states 

Proposition VI. 1. The states lu„ are locally quasiequivalent to a Hadmard 
state if n > N . Furthermore, the difference between the twopoint functions fJn 
of LOn and those of a Hadamard state is given by a C^-^+^ kernel. 

Proof. According to Thm. IVI.ll an adiabatic state of infinite order (defined, as 
described above, by a symbol 6 e S (/,R)) is Hadamard, so it is sufficient to 
show that adiabatic states of order n (described by a symbol b e S"(/,]R)) are 
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locally quasiequivalent to such a state. Now b ~ bn is by definition a symbol of 
order — (n + 1), so in particular, 

\\bks - bnksW < c|l + k\-^-' for all fc > 0. 

The criterion on local quasiequi valence, Thni. IIII.ll then immediately proves 
that the states are locally quasiequivalent ii n > N. 

Let An € OP"""""^ (S(t), I?M) be the operator associated to the symbol 
a„ = b — bn via Eq. {T)) at some t £ I. By standard theorems, e.g. in PUI II, 
Prop. 2.7], the associated kernel on S(t) x E(i) is in (7"-^+^ for n > iV - 1. 
The difference of the twopoint function of an adiabatic state of infinite order 
and one of order n is 

^W(/,M-^i±)(/,M = ±((^/i) rs(i),A.(^/) rs(i)>f 

Now the causal propagator ^ propagates m times differentiable initial data to m 
times differentiable solutions. Therefore the above difference must (n — iV + 1) 
times differentiable in M x M. D 

Theorem VI. 1. The adiabatic states are Hadamard in the sense of Def. lIV.H 

Proof. In order to prove that f?^^-* has the wave front set described in Def. lIV.fl 
we shall employ the following result due to W. Junker TF, Thm 3.12] . This result 
has originally been obtained for scalar fields, but a careful analysis of the proof 
shows that it can be adapted to the spinor case. We present here a modified 
version which is tailored to our situation. E is the causal propagator for the 
spinorial Klein-Gordon operator P. 

Theorem VI. 2. Let Q{t) he an elliptic PDO on S(t). Let I be an interval 
containing t and A± G OP(I] x L,DM) such that there exist PDO's R± £ 
OP(S X /, DM) which have the property R±{in'^'V ^ + A±) — P modulo smooth 
and 

Qb± c{(x,C)er*Af\{0}ie>(<)0}, 

where Q_r^ is defined in Eq. (|38|l . Then the spinorial bidistributions 
fS^\h, f) = ((zn^V^ + A±{t))Eh, Q{t){in^W^ + A±{t))Ef)t 

have wave front set WF'(^(='=)) C C^^^ 

We apply the lemma to A^ = —T and A+ = T* and Q as in the definition 
of the twopoint functions, Eq. (|77|l . Then Eq. (j^ . and Eq. (J5HJ| provide us 
with operators 

i?_ = -(m'^V,, + iK + T), R+ = -{inf'XI ^ + iK -T*) 

as in the statement of the above theorem. Clearly, since cri(T) — ai{T*) — 



Qr^ = K^, I n^i^ = ±V-h'"'UA = {{^, I e > (<)0} 

Noting that {—if7 — m)E = ^ and using the fact that the wave front set cannot 
become larger upon acting with a PDO on a distribution, we can apply Junker's 
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theorem to G^^'> (given by Eq. ^) and obtain WF{f!^^'>) C C^^K It remains 
to show equahty in the above inclusions. The anticommutation relations imply 
that ^(+' + ^^^^ — i^. If the causal propagator ^ had wave front set W — 
{(a;i,Ci,a;2,6) I (a;i,^i) -- (2^2,6)} = C^^^ U C^-^ (and this will indeed be 
shown) then 

W = WF'(^) c WF'(^W) U WF'(^(-)) c C(+) U C(-) = W, 

thus in fact equality would hold in the above inclusions. We have to show that 
jS" has indeed wave front set W. From ^'^ = ^ and the antilinearity of the charge 
conjugation it follows that 

WF'(^) = -WF'(^). (35) 

One has (see 0]), jS* f S x E = ifil where 1 means the identity on the Cauchy 
surface. From T^, Thm. 2.22], one knows that 

WF'(^rSxS) C d<^*oWF'(^):-{(xi,rf0^^(a),X2,d4.(6))l 
((/.(xi), a, 0(2^2), 6) e WF'(^), x,,X2 e E} 

where (p : T, ^ M is the embedding map. Now assume that (3:1,^1,2:2,^2) is in 
W but not in the wave front set of ^. By the propagation of singularities, we 
can assume that there is an element (a;, £,,x,S.), x € E, which is in W but not in 
WF'(^). By Eq. (ESJ, also (x, -^,x, -^ ^ WF'(^). Since ^ must be a nonzero 
null covector, it is impossible that the nonzero element {x,d(j)l.{^),x,d(j)l.{^)) is 
in d(t)* o WF'(jS') D WF'(i/il). But the latter set is actually equal to 

WF'(5W) = {(x,e,2;,OieeT*E\{0}}, 

and so must contain any element of that form, a contradiction. D 

VII Comments 

1. At adiabatic order zero we get Bo(t) = ^\H{t)\-'^{H{t) + \H{t)\). i.e. 
Boit) projects on the instantaneous positive frequency solutions at time 
t. Clearly, if R is constant in a neighbourhood of i, then there are no 
further corrections to this operator at higher adiabatic orders, and Boit) 
defines a pure quasifree Hadamard state. If R is not constant near t, then 
the operators Bn{t) give states which are not Hadamard in general but 
only reproduce the highest order singularities of the Hadamard form, in 
the sense that^ (modulo C"-^+^), 



fJ<t\xuX2)^ 



± {if + m) 



{N-l)/2 n~N+2 



Uj{xi,X2)cr^l + 2_] Vj{xi,X2)cr^loga±e 



Hence, for n > N , adiabatic states will allow for a point-splitting renor- 
malisation of the stress-energy tensor T^^, as described e.g. in 34 (the 



^Note that the next term in this series would be a'^ N+3 log cr, which is (n — N + 2) times 
differentiable. 
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difference to a Hadamard state must be at least in C^, because the stress 
tensor contains 1 derivative), but such of lower order will not in general. 
In other words, we see that any adiabatic vacuum state which allows for 
a point splitting renormalisation of the stress-energy tensor will be locally 
quasiequivalent to a Hadamard state. 

2. In 121], the authors observe that the expected stress tensor of a Dirac field 
diverges for the 'energy-minimising states' proposed in that work. This 
is explained by our analysis, since their states are simply the adiabatic 
states of order zero just discussed in diguise. 

3. Following the strategy of Liiders and Roberts [22] for the scalar Klein- 
Gordon field, 23 proposes another definition of adiabatic states for the 
Dirac field. We do have some doubts as to whether their definition really 
yields a positive state, moreover the role of positive and negative frequency 
modes remains obscure in 135^. It is therefore difficult to see how their 
definition relates to ours. In view of the analysis carried out in this paper, 
we do not believe that the states proposed in [SS] are of Hadamard type, 
as suggested by the author. 

4. It is possible to construct a Hadamard state in a general globally hyper- 
bolic spacetimes along the same lines as in the previous section if one can 
construct a hermitian, positive PDO Q such that 

L+QL*^+L*_QL^ = 1. 

However, unlike in the simple case of a RW-spacetime, the construction of 
such a Q seems to be harder. We are currently working on this problem. 

VIII Appendix 

VIII. 1 Spinors on flat space and representation theory 

In this appendix we find the (generalised) eigenfunctions of the spatial Dirac 
operator ^ on M.^ . To this end, we first write this operator in polar coordinates. 



^^ = {^e + ^^) l''^ + l^N-ii^^ 



where ^ ]\j_i is the Dirac operator on §^ ^ . The eigenfunctions of this operator 

EH, 

Vn-1^1^^ = Ml + iN- l)/ml^J, ; = 0, 1, . . . , m = 0, 1, . . . , dz, 

may be used to find the spectral decomposition of Dirac operator on R^ . We 
first set 

£(±) - 1. ff (-) ± ,yNA+) 

In order to find the eigenfunctions of )(7 we insert the ansatz 

Xm™± (0, ^) = c(kl) [au mt} (f^) ± *&fci (^)eL^ (^)) , f^ e S^-i (36) 
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into 

WXklms — ~k Xklms- 

This leads to the differential equation 

'2 N-l^ l{l + N-2) ,2' 
dn + — T^de + T7^ + k^ 



akiiO) = 0, 



for aki (and similarly 6^;). The unique regular solutions to these equations are 
given by Bessel functions, 

bkiie) = 0-(^-^^/'Ji+N/2ike). (37) 

The normalisation factor in Eq. (|36(l is determined from the condition 

{XklrnsiXk'l'm's') = S{k — k )Sll' SmTn' Sss' , 



one finds c{kl) — y^k/2. In this work we also need the spectral function 
(Plancherel measure) defined by 

^^(fc) = E^«-(o)^xi:L(o). 

Im 

From the expression Eq. 1)37(1 and behaviour of Bessel functions at 6* = it is 
seen that only the term with / — will contribute, leading to the result 

PN{k) 



2^/2vol(S^-i)r(A^/2)2' 

VIII. 2 Notions and results from microlocal analysis 

For convenience we mention some results and definitions from the theory of 
distributions and the theory of pseudodifferential operators (PDO's). If not 
indicated otherwise, these may be found in standard textbooks, for example 
|3()| ll3j. PDO's generalise ordinary differential operators in the sense that they 
give meaning to fractional powers of derivatives. They are defined in terms of 
so-called symbols. We shall not give the most general definition of a symbol 
here, since only a certain class of symbols is important for this work. 

Definition VIII. 1. Let ^ be a subset of E" and to be a real number. Then 
a symbol of order to, is a function a G C°°(£?,R") such that for every compact 
subset K oi ff the following estimate holds 



ni-|/3| 



D^Dla{x,0 <aA/^(l + |^|) 

for all multiindices a,p. D" is i'^'c?"^ . . .9"". The set of all such symbols is 
denoted by S'"(^,]R") and one also writes S"°° = flmS". 

There is the notion of the asymptotic expansion of a symbol which is an 
important tool for constructing PDO's. Suppose Gj e S'"^(i^,R") for j = 
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0, 1, 2, . . . with nij monotonously decreasing to minus infinity. Then there exists 
a e S™''(^,R") such that for aU N 

N 

i=o 

and a is defined modulo S"°°. One writes a ^ J2j «j- If a G S'"(^,R") then 
the operator 



Au{x)= / e"^a{x,^)u{0 



(27r)" 

is said to belong to 0P™(€?), the PDO's of order m. A is a continuous linear 
operator from &{ff) to C°°(R"). By the Schwartz kernel theorem it is thus 
given by a distribution kernel Ka E i^'(^ x ff). Ka is smooth off the diagonal 
in ^ X ^ and smooth everywhere in ^ x ff ii A E 0P^°°(^). Hence the asymp- 
totic expansion of a symbol uniquely determines a PDO modulo smoothing 
operators. The above statement carries over to matrix valued symbols without 
major changes. A principal symbol am{A) of A e 0P™(^) is a representor of 
its symbol in S'"(^)/S'"~^(^). It can be chosen such that it transforms con- 
travariantly under a change of coordinates (giving thus a well-defined function 
on the cotangent bundle) and it behaves multiplicatively under multiplication of 
two PDO's. On a manifold M (or more generally on a vector-bundle E) PDO's 
are defined to be the continuous operators on ^(M, E) which have the above 
properties in each coordinte patch. 

We come to the definition of the polarisation set of a vector-valued distri- 
bution u = (u\ ■ ■■ ,u'') € &'{ff)'', ff an open subset of K". For details of the 
definition and the subsequent results see |3]- 

Definition VIII. 2. The "polarisation set" WFpoz(w) of a vector- valued distri- 
bution u is defined as 

wFpo,(")= n ^^' 

Aeov°, Auec°° 
where 

^ = {{x,tw) e T*e X C*^ I aQ{A){x,i)w = 0}. 

From the transformation properties of the principal symbol it is clear that the 
definition can be carried over to the case of distributions with values in a vector- 
bundle E. WFpo/(u) is then seen to be a linear subset of ■k*E, tt : T*M -^ M 
being the canonical projection in the fibres of the cotangent bundle. The "wave 
front set" WF(m) of a distribution is obtained by taking all points (x, S,) G T*M 
such that the fibre over this point in WFpo;(u) is nontrivial. The microlocal 
properties of the bidistributions considered in this work are more conveniently 
described in terms of their primed polarisation set, WF ;, which is obtained 
from the usual one by reversing the sign of the covectors in the second slot (the 
primed wave front set is defined similarly). 

There is an important theorem on the polarisation set of distributions u 
satisfying Pu £ C°" for differential operators P of real principal type, which 
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goes under the name 'propagation of singularities' |31E|- Such operators are 
defined as follows (in the following we discuss the simple case where P acts in 
the bundle E — M x C'', but all results can be generalised in the obvious way 
to nontrivial vector bundles): 

Definition VIII. 3. A k x k system P of differential operators on a manifold 
M with principal symbol poix, ^) is said to be of real principal type at (y, r/) if 
there exists a k x k symbol po{x,^) such that 

Po{x,£,)Poix,£.) = q{x,£,)^k 

in a neighbourhood of (y,??), where q{x,S.) is scalar and of scalar real principal 
type, i.e., d^q{x, £,) ^ for all ^ 7^ 0. 

One sets 

Qp = {ix,O\detpoix,O = 0}. (38) 

If / is a C°° function on Qp with values in C''', then one defines 

Dpf = X,f + ^{po,Po}f + tpop'f, (39) 

Xq being the Hamiltonian vector field of q, 

Xq = dxqd^-d^qdx, {pi3,Pc} = d^Po d^Po - dxPo d^po , 
p'' = pi + -d^Dj;PQ, <j{P) -^ Pq + Pi + P2 + 

One can prove that Dp is a partial connection along the Hamiltonian vector field 
restricted to Qp. Since there is some arbitrariness in the choice of the symbol 
p, the partial connection is not uniquely defined. One can however prove that 
the remaining arbitrariness is irrelevant in what follows. 

Definition VIII. 4. A Hamilton orbit of a system P of real principal type is 
a line bundle J^p C ^p \ c, (c is an integral curve of the Hamiltonian field on 
Qp, c{t) = Xq{c{t))) which is spanned by a sections / satisfying Dpf = 0, i.e. 
Jffp is parallel with respect to the partial connection. 

Tiieorem VIII. 1. Let P be of real principal type and u a vector-valued distri- 
bution. Suppose {x,S,) ^ WF(Pm). Then, over a neighbourhood of{x,^) in Qp, 
WFpoi{u) is a union of Hamilton orbits of P. 
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